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Abstract 



It is well known that the electroweak phase transition (EWPhT) in extensions of the 
Standard Model with one real scalar singlet can be first-order for realistic values of 
tH- \ the Higgs mass. We revisit this scenario with the most general renormalizable scalar 

■^ ■ potential systematically identifying all regions in parameter space that develop, due to 

tree- level dynamics, a potential barrier at the critical temperature that is strong enough 
t^ \ to avoid sphaleron wash-out of the baryon asymmetry. Such strong EWPhTs allow for 

a simple mean-field approximation and an analytic treatment of the free-energy that 
leads to very good theoretical control and understanding of the different mechanisms 
that can make the transition strong. We identify a new realization of such mechanism, 
IT/ \ based on a flat direction developing at the critical temperature, which could operate in 

'j_j ■ other models. Finally, we discuss in detail some special cases of the model performing 

. 5t , a numerical calculation of the one-loop free-energy that improves over the mean-field 

approximation and confirms the analytical expectations. 



1 Introduction 

The search for physics beyond the Standard Model (SM) has strong theoretical and experi- 
mental motivation and the simplest extension is to enhance the SM by a scalar gauge singlet 
degree of freedom. This minimalistic model (and its cousins with a complex singlet or su- 
persymmetric versions of it) can be very successful in explaining various phenomena that 
cannot be explained by the SM: dark matter [I]- [8], spontaneous B — L breaking [9]-[Ti] 
and the baryon asymmetry of the Universe [T5| [16], often leading to characteristic collider 
signatures [T7]-p3]. 

One prominent difference between the SM and its singlet extensions is the following. 
While in the SM the LEP bound on the Higgs mass {Mh > 114 A GeV [23]) implies that the 
electroweak phase transition (EWPhT) is not first-order but a smooth crossover j25j, the 
addition of a singlet can lead to strongly first-order EWPhTs [26]-|l2] for realistic values of 
the scalar masses. Moreover, with such strong EWPhTs, not only the observed baryon asym- 
metry of the Universe can be explained through electroweak (EW) baryogenesis (provided 
the model also contains additional sources of CP violation) but the process of EW symmetry 
breaking can also leave the trace of a stochastic background of gravitational waves [43j. 

The aim of the present work is to revisit the EWPhT in the most general renormalizable 
extension of the SM with one additional real scalar singlet. Although this issue has been 
studied in the past [26]- [12], (both numerically and analytically at different levels of general- 
ity), we believe that a thorough analytical understanding of the rich spectrum of possibilities 
for a strong EWPhT this model offers is still lacking in the literature. The analysis that 
comes closest to this task is Ref . [35] , over which we will improve in a number of aspects. 

In the SM and some extensions of it, a first-order EWPhT is caused by the thermal effects 
of bosons coupled to the Higgs, that generate a cubic term in the Higgs scalar potential. 
Although this can be successful in many cases, it requires sizeable couplings of these bosons to 
the Higgs and the effect can be screened by thermal masses when daisy resummation is taken 
into account. In this article, we concentrate on EWPhTs for which the barrier separating 
broken and symmetric vacua is not generated by the previous thermal cubic correction but 
rather by tree-level effects. These tree-level effects lead in general to stronger EWPhTs as 
Vc = v(Tc), the Higgs vacuum expectation value (VEV) at the critical temperature Tc (that 
controls through the celebrated ratio Vc/Tc the sphaleron erasure of the baryon asymmetry). 



is now proportional to some T-independent dimension 
Vc/Tc can become potentially very large for small T^ 







ul parameter in the potential; hence 
The parameter space of this model 



^Strong EWPhTs are particularly welcome if, as suggested by [55], magnetic fields generated at the 



is quite rich and we will see that these tree-level barriers are not necessarily related to the 
presence of cubic terms in the potential, as is often assumed. 

We begin in Section |2] by studying the structure of the tree-level scalar potential of the 
model. In particular we are interested in potentials where the EW breaking and preserving 
minima are degenerate, since this is the situation that arises for strong phase transitions. 
Indeed, in this case it is well justified to use the mean- field approximation for the free-energy, 
which will have the same structure as the tree-level potential (with temperature-dependent 
parameters). Differently from previous analyses, we will introduce a novel set of parameters 
particularly convenient for the discussion of the vacuum structure of this kind of potentials 
(and which might also be of use for phenomenological studies of the scalar sector of these 
models). In spite of the large number of free parameters (eight) we have to deal with, this 
new parametrization will allow us to identify very easily, and analytically, the structure of the 
potential: its stable minima and the existence of a barrier between them. As a result, beside 
developing a better understanding of the ingredients necessary to get a strong EWPhT in 
this model, we will find new scenarios with strong EWPhTs that had not been identified 
before (involving in particular flat directions at the critical temperature). 

In Section |3] we discuss thermal corrections to the scalar potential and explain our strategy 
to search for regions in parameter space with strong EWPhTs, which we summarize in 
Table [1] The idea is to start from a potential with degenerate broken and unbroken minima 
and a barrier between them: this is identified with a potential at Tc which gives a strong 
EWPhT. Its parameters are then evolved to lower T to find their values at T = 0, where 
they can be used for phenomenological purposes. After identifying the regions in parameter 
space that give a strong EWPhT, we then perform a more precise analysis including one- loop 
corrections (at T = and finite temperature) without resorting to high-T expansions and 
including daisy resummation. Although the strength of the EWPhT is somewhat reduced 
with respect to tree-level, one still gets sizable values for it. Our results confirm then the 
expectations based on the tree-level analysis. 

In the rest of the paper we apply these tools to particular realizations of the model 
previously considered in the literature: the Z2-symmetric case in Section HJ a particular 
supersymmetric incarnation in Section |5l and a case with a very light scalar in Section [61 
Finally, we study some examples of the general case in Section [7] and conclude in Section |H1 
Appendix |A] contains some technical details of the full one- loop analysis, including the T = 
renormalization conditions. 



EWPhT lower the sphaleron energy so that larger values of Vc/T^ are required to avoid baryon washout. 



2 Tree- level Scalar Potential 

2.1 Parametrization of the Potential 

To begin with, we look for a convenient parametrization of the potential that ensures control 
over its minima: this will allow us to easily identify which are the global minima and whether 
or not they are stable. The most general (renormalizable) tree-level potential for the SM 
Higgs field h and the singlet s depends on 8 parameters, 

V = --^4^'' + -Xhh" + -^Jiy + -Ks^ + -^lmsh'' + -XmS^h'' + ^i\s + -/igs^ . (2.1) 

Note that a redefinition of the singlet field by a constant shift, s — )■ s + cr, simply amounts 
to a redefinition of the parameters fif, ^1, fil, ^3 and fim but does not change the physics, 
being just a coordinate change: 

nl -> fil + Asff^ + fi?y + Ilia , 

/^s -^ f4 + SXgCr^ + 2/i3(T , 

/W3 -^ /Lt3 + 3AsCr , 
/^m -^ /^m + 2Am(J . (2.2) 

This shift is often used to get rid of one of the initial parameters, choosing e.g. /ii = or 
A^s = Oo However, different choices can be useful in different circumstances so we refrain 
from any particular choice at this stage and keep the discussion as general as possible. Still, 
it is of advantage to choose a parametrization in which the shift symmetry is realized in a 
more explicit way. Beside this property, the new parameters will allow for a more direct 
theoretical control of the structure of the potential. The parameters we introduce are: the 
vacuum expectation values v = (h) and w = (s) in the broken minimum; the three elements 
of the scalar squared-mass matrix, evaluated at the broken minimum (as indicated by the 
subscript b), 

(2.3) 
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the mixed quartic coupling A.^,; the effective coupling A^ defined by 



A' = A,A, - ^Xl , (2.4) 



^If the potential is invariant under the discrete Z2 symmetry s -^ —s then /i„i = /^i = /ia = is the best 
"coordinate frame" for the singlet as it makes explicitly apparent such symmetry. 



that appears recurrently in different contexts (note that A^ can be negative); and finally, the 

combination 

m* = X'^w + -Xhf^s - ^Xmfim , (2.5) 

o o 

which can be checked to be shift-invariant. 
To sum up, our parameters are 

{v, w, ml, ml, m],^, A^, A^ mj . (2.6) 

With the exception of w (transforming as w — )■ w — a), all these parameters are shift- 
independent. For reasons that will become clear when we discuss the thermal potential, it 
is convenient not to assume at this point that v takes its standard value vew = 246 GeV. 
To avoid confusion later on, we reserve the notation vew for the latter value. The relations 
between the old and the new parameters are: 

^^h = Ij^h + -"l^sh - ^^mW , (2.7) 

^il = ml- -XmV^ + ^-2 [-2Xmmlf^ - 8m,v + {AX^ + Xl)vw] , (2.8) 

3v 
/^3 = ^^ [4^m,v + Xmml^ - {AX^ + Xl,)vw] , (2.9) 

9 

llm = 2^^~2XmW, (2.10) 

V 

2 

1^1 = -{XmVW-ml^)-mlw + ^-^[3X„,ml,^ + 12m,v-iAX^ + Xl)vw] , (2.11) 
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A. = i^, (2.12) 



.2 , ^2^ V"" 



Xs = (AX' + Xi)^^. (2.13) 

This change of variables is non-singular: indeed its Jacobian is simply 12{mlm\ — m'^f^)/m,f^ 
which, as shown in the next Section [eq. (12.151) ]. is always positive and non-singular. 

The potential in the new parametrization reads (up to an appropriate constant) 

V = <{e-v'y + ^{h^-v'){s~w) + \[2ml + XUh'-v')]{s-wr 

+ —^{Xmmlf^ + Am,v)vis-wY + —^iAX^ + Xiyis-w)\ (2.14) 

2mi 8mi 

where we have expressed it as a polynomial in {s — w), showing explicitly that now the 



parameter w can be used to absorb any shift in s, leaving the other parameters of the 
potential invariant. 



2.2 Structure of the Potential 

As a good starting point for the analysis, we would like our tree-level potential to be well- 
behaved. First, it should have a stable broken-phase minimum {v,w). This is guaranteed by 
using V and w as input parameters and by a judicious choice of the mass-matrix elements in 
02 .Sp . such that both mass eigenvalues are real and positive, that is, 

Det Ml = mlml = mlml - m^^ > . (2.15) 

Alternatively, one can use directly the mass eigenvalues, m^^ , m1^ and the scalar mixing 
angle, ash, as input parameters, imposing the relevant experimental bounds [44] (suitably 
modified to take the singlet mixing into account). 

2.2.1 Stability 

We also want that the tree-level potential does not have unbounded-from-below directions. 
The large- field behaviour of the potential is dominated by the quartic part, which in our 
parametrization reads: 

^4 = ^^ [{mlh'' + Xmv^s^f + 4AVs^] . (2.16) 



For Am < 0, the squared-sum term vanishes along the directions h = ±y/—XmVs/mh, so 
that A^ > is required to ensure stability. For A^ > 0, on the other hand, the squared-sum 
term is always positive and the positivity condition on A^ can be relaxed but, to maintain 
stability along the s-direction, one should require A^ > — A^/4 in that case [so that A^ > 
in the parametrization of eq. (12. ip . 

The fact that {v,w) is a local minimum does not guarantee that it is the global one: a 
deeper minimum might exist. We will derive below the necessary and sufficient conditions 
for this situation to arise. These conditions turn out to be extremely simple in terms of our 
new parameters. 

2.2.2 Local Minima at h ^ 0. 

Now we search for further local minima of the potential in order to ensure that {v, w) is the 
global one or, at least, degenerate with an unbroken minimum (as relevant for the EWPhT). 
The stationary points of the tree-level potential lie on the curves along which dV/dh = 0, 

^ = ^ |/i = 0, and h^ = Dl{s) = ^^{2fil-fL^s-X^s^)\ . (2.17) 
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Figure 1: The curve /i^ = -D^(s) in the (/i^/f ^, s/w)-plane. The different cases correspond 
to: Am = (upper left); Am > (upper right); Am < with |Am| < "^\hl K"^hP'^^ (lower left); 
Am < with |Am| > ^h/ (''^1'^^) (lower right). The unphysical region h^ < is shaded gray. 
The EW breaking minimum at h"^ = v"^ and s = w is marked by a black dot. 



In our parametrization, the curve Df^{s) reads 
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(2.18) 



In particular, we have D\{w) = v"^ (by construction) and it is interesting that this curve is 
independent of the parameters m*, ml and A^. Its shape, determined by Am and mil^, will 
be relevant later on. Fig. [1] shows the different possibilities for the -D^(s) and h = lines 
in the {h'^/v'^, s/w)-plane, with a minimum at s = w as indicated. These lines separate the 
plane in regions of definite sign of dV/dh. 

For Am = 0, the curve -D^(s) is a straight line, intersecting the axis h = at one single 
point (Fig. [H upper left). For the special case in which also ?7i^^ = one simply has 
Dl{s) = f^, and the corresponding line is parallel to h = 0. 



In general, with Am 7^ 0, -D^(s) is a parabola and can have three qualitatively different 



forms. If Am > 0, -D^(s) curves down and intersects h = at two points. If Am < 0, -D^(s) 
curves up. For |Am| > "^sh/l"^^"^^)' ^hi^) ^'^^^ ^^^ intersect h = while, in the opposite 
case, it has two intersection points. For either sign of Am, the two intersection points with 
the s-axis are given by 

1 



s± — w 



-^Ih ± V "^sh + ^mv'^ml 



(2.19) 



AmW 

All these different possibilities are illustrated in Fig. [1] We can use the shift freedom to move 
the axis of the parabola to s = 0, which corresponds to enforcing 

2 

w = p^, (2.20) 

Amf 

or to setting ^^n = in the original parametrization. This is the choice we generally adopt 
in our plots. 

Next, we consider possible additional stationary points along the curve h"^ = D\{s). The 
potential along such curve, V[Dh{s), s], is a quartic potential in s that can be minimized in 
a straightforward manner. More explicitly, the minimization equation dV[Dh{s), s\/ds = 0, 
leads in the general case to a cubic equation of the form 

a{s - wf + b{s -wf + c{s-w) + d = 0, (2.21) 

with 

a = 2X\'^, b = 6my, c = DetMl, d = . (2.22) 

The nature and number of real solutions this cubic equation has is determined, as usual, by 
the discriminant 

A = ISabcd - 4b^d + b'^c^ - 4ac^ - 270^^^ . (2.23) 

For A < there is only one real root, corresponding to a single minimum, the electroweak 
one; for A > there are three real roots (the previous minimum and two other stationary 
points); for A = the two additional roots merge in an inflection point. Notice that these 
additional stationary points are only physically relevant if they appear in the region with 
Dl{s) > (the interval [s_,s+], with s± defined by eq. (I2.19p . if Am > 0; or the intervals 
[— oo, s_], [s+, oo] if Am < 0). With our coordinates, 

sign(A) = sign [9mlv^ - 2A2DetA^^] . (2.24) 

Recall that Det A^f = mlml—m^f^ > from (I2.15p . This means that the necessary condition 
to have an additional stationary point along -D^(s) is 

7 



where A^ > 0. 

In the case < A^ < A^, the two additional stationary points are another minimum and 
a maximum separating it from the EW breaking one. Their location is also easy to obtain: 
they appear at {Dh{w±),w±) with 



3m 

w± — w = 



± 7rr^^/9mlv^ - 2A2DetM2 . (2.26) 



2A2 2\-^v 

By evaluating the potential at these points, it is straightforward to obtain the condition for 
the minimum at (f , w) and the additional one from eq. fl2.26p to be degeneratejj 

A2 = 8AV9 . (2.27) 

For A^ < [which requires A^ > from the stability discussion below eq. fl2.16p ]. the two 
additional stationary points are two maxima, with the EW minimum between themcl. The 
EW minimum will still be the deepest (physically relevant) minimum along -D^(s) provided 
V{v,w) < V{0,s±). If this is not the case, a deeper minimum must exist at /i = 0. We 
discuss such situation in the next Subsection. 

To summarize, as illustrated by Fig. [21 for a potential with all parameters fixed except 
A^, the possible stationary points away from h = Olie along a fixed curve -D^(s), independent 
of A^, with a minimum located at {v,w) by construction. For large enough A^ > A^ > 0, the 
minimum at {v,w) is the only stationary point. When A^ = A^ an inflection point develops, 
while for smaller A^ there are two minima. For A^ = 8A^/9, the new non-standard minimum 
is degenerate with the one at {v,w) and for smaller A^ our minimum {v,w) is no longer 
the lowest one. Whether these other minima are physically relevant or not will depend on 
whether they appear at positive values of -D^(s) or not. 

2.2.3 Local Minima at h. = 0. 

To discuss the possible presence of minima in the direction h = 0, which might be deeper 
than the electroweak vacuum, two simple facts are relevant. First, it is useful to note the 
following relation: 

2 

m 



V[D,{s), s] - V[0, s] = --iDlis) < . (2.28) 



^y2 



•^In this degenerate case, w- — w = 2{w+ — w), corresponding to a potential V[Dh{s), s] symmetric under 
(s — w+) — >■ — (s — w+) (this is generic for a quartic potential with two degenerate minima). 

^In this case V[Dh{s), s] is unbounded from below for s — >■ ±oo; however, this is not a problem because 
this region is not physical when A„i > 0, as discussed above eq. (|2.24p . 
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Figure 2: Example for the dependence of the potential on A^. The upper left plot shows 
Dl{s) (solid blue) and D'^{s) [along which dV/ds = 0, see eq. (12.351) ] (solid red) for several 
values of A^: A^, 8A^/9 and A^. The intersections of these two curves correspond to the 
stationary points of the potential. The remaining plots show the potential along Dl{s), 
D^{s) (same color coding) and h = (dashed) at the indicated values of A^. 



A glance at Fig. [T] explains this fact as a result of the sign of dV/dh in the region between 
the lines h = and /i^ 



Dl{s). Second, 
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(2.29) 



This implies that minima along h = can only appear in regions with -D^(s) < Oo Moreover, 
in order to locate such minima, it is enough to minimize l^[0, s]: if a minimum found in this 
way appears at -D^(s) < 0, then (I2.29P ensures that it is also a minimum along the h- 
direction. We also conclude that, for fixed s the potential can have only one minimum: at 
h = Oior Dl{s) < or at /i 7^ for Dl{s) > 0. 

More can be said about these minima at /i = by considering the following general 



^An alternative way to see this is to note that d^V/{dh)^ = along the line h^ = Df^{s)/3, which divides 
the (h^, s)-plane in two regions with opposite signs of d^V/{dh)^. 



statement about the tree-level potential: // there are two local minima with h"^ > then 
there is no local minimum with h = 0. To prove this, assume there are two local minima 
with h^ > and one with h = 0, chosen to be the absolute minimum along h = and to lie 
at s = 0. Consider a curve of the form 

h"^ = D^{s) = as + /3s'^ = (3s{s-si) , (2.30) 

that passes through all three minima (two minima with equal s are not possible and so, 
this curve always exists). The potential along this curve, Vd(s) = V[D{s),s], is a quartic 
polynomial and, therefore, can have at most three stationary points. Two of them are the two 
minima with /i^ > while the minimum at /i = is in general not a stationary point of Vd{s). 
The reason is that dV[D{s), s]/ds = [dV/d{h'^)]dD^{s)/ds + dV/ds, but dV/dh = does not 
in general imply dV/d{h'^) = because /i = leads to dV/dh = even if dV/d{h'^) ^ OO 
Furthermore, minima must be separated by maxima and in most geometrical arrangements 
of the locations of the three minima this will require more stationary points than the allowed 
maximum of three. The only non-trivial case occurs if the minimum at /i = is between the 
other two (ordered by their s-coordinates) and si 7^ 0, but in that case it is straightforward to 
see that Vd(si) < Vd(0), contradicting our assumption that s = is the absolute minimum 
along h = 0, and this concludes the proof O 

On the other hand, if there is only one local minimum with h'^ > there can be up to two 
local minima with h = 0. The case with two minima at /i = requires that the potential 
Vd{s) has one minimum (the one corresponding to /i 7^ 0) and is negative for s — )■ ±00. 
In addition, the two minima at /i = lie at both sides of the EW minimum. (An explicit 
example will be given in the Z2 symmetric case below). 

Having these facts in mind, we can compare minima a.t h = with our minimum {v,w) 
and discuss what are the conditions on the parameters of the potential for {v, w) to be the 
global minimum. We will illustrate this with an example in Fig. [2l which plots the potential 
along h"^ = Dl{s) and h = for different values of A^. Consider first the case in which {v, w) 
is the deepest minimum of V[Dh{s),s] (i.e. we have A^ > 8A^/9), then ( 12.28^ immediately 
implies it should also be the global minimum of V[h, s]. Such case is shown in Fig. [21 upper 
right plot. Cases with A^ < 8A^/9, for which a deeper minimum along -D^(s) appears (or, 
if A^ < 0, when V[Dh{s), s] is unbounded from below), might still have {v,w) as the global 
minimum (this happens if the new minimum is in the unphysical region, Dl{s) < 0), as in 
Fig. [21 lower left plot. In such cases one needs to check the minima along h = 0, which might 



^A relevant exception to this is h^ — Df^{s) as then we always have dV/d{h'^) = by construction. 
^Obviously, this "theorem" applies to the tree-level potential only and can be violated through loop 
corrections. It is nevertheless useful in order to identify large barriers created by tree-level effects. 
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be deeper than the EW one without violating fl2.28p . As A^ gets more and more negative, 
minima along h = might become the global minimum. In general, when all parameters 
except A^ are fixed, the potential along h = decreases with decreasing A^ [see the explicit 
potential in eq. f l2.14p ]. Hence, there is a definite value A^ for A^ that separates the region in 
parameter space in which the EW minimum is the deepest one from the one where it is not. 
It is clear that — A^/4 < A^ < 8A^/9. The exact value of A^ will be determined below by 
requiring degeneracy between the EW minimum and the second minimum [see eq. fl2.43p ]. 
an example of which is shown in the lower right plot of Fig. [21 Studying such degenerate 
cases will be very relevant for the discussion of strong phase transitions in the next Section, 
so we turn to this issue next. 

2.3 Coexisting and Degenerate Minima. 

The most interesting cases for the phase transition study are potentials with two degenerate 
minima: the EW-breaking one, at {v,w), and the symmetric one at {0,wq). While we could 
use the shift of eq. (12. 2 p to specify the value of the singlet field VEV in the unbroken phase, 
''^o {wq = is often used in the literature), here we will keep again the shift-invariance 
explicit since, in our parametrization, such choice would simplify intermediate expressions 
only marginally. 

We will next show that, out of the eight initial parameters, only three have an impact 
on the shape features of the potentials with two degenerate minima. Two parameters can 
be removed thanks to the shift-symmetry and the requirement of degeneracy of the minima. 
Two more parameters can be removed by rescaling the potential in the s and h directions. 
Finally the overall scale of the potential has no qualitative meaning in this discussion, leaving 
us with three parameters. In the following we present a parameter choice, a refinement with 
respect to eqs. ( I2.7p - (]2.13p . which is especially handy in describing the qualitative features 
of a potential with degenerate minima. 

To discuss the minima, let us consider the curves dV/dh = and dV/ds = more 
systematically. The curves h = and Dl{s), at which dV/dh = 0, were already introduced 
in eqs. fl2.17p . We begin with A^ 7^ 0, in which case we can rewrite the curve -D^(s) as: 

DUs) = ¥-^{s-w,)\ (2.31) 

mi 



11 



where, using both the notation of eqs. (I2.7p - fl2.13p and the original notation, 

2 

"^v = ^ - v^ = -vr ■ ^^-^^^ 

In the plots, we will generally choose our singlet coordinates to have Wp = (setting /x^ = 0), 
so that the parabola -D^(s) has its axis at s = 0. As discussed in Section |2l the minimum at 
h = can only be located in the regions with DKwq) < 0. 

The potential is a quartic in s and hence has at most three extrema with respect to s for 
fixed h (out of which at most two are minima). For general values of h, 

OT 7" -1 

— = -h\fXm + 2\mS) + ifxl + /i^S + /XgS^ + A,s^), (2.34) 

and the curve dV/ds = in the (/i^, s)-plane as a function of s is given by 

e = Dlis) ^ _^d±fk±J^^f±M , (2.35) 

or, in our parametrization: 

( g — tv] 
Dl{s) = v^ - % _ — r [2mlml + 3i;(A^m^^ + 4m,t;)(s - w) + {AX'^ + A^)(s - wf] . 

(2.36) 
This function is single-valued for fixed s and has a pole at s = Wp. When this pole is 
canceled by a zero of the numerator, the line s = Wp is also a solution of dV/ds = (this 
is e.g. the case for a potential with a Z2 symmetry s — )■ —s) and enters the discussion. The 
asymptotic behavior of -Df (s) at large s is D'^{s) — t- — 2AsS^/Am and hence is qualitatively 
different depending on the sign of A^. In the following, we distinguish four different cases 
given bjo 

(a) A„ < , (wq- Wp){w -Wp) > , (2.37) 

(6) A„ < , {wo- Wp){w -Wp) <0 , (2.38) 

(c) A„ > , {wq- Wp){w -Wp) <0 , (2.39) 

(d) Xm>0 , {wo- Wp){w -Wp)>0 . (2.40) 

The sign of A^ determines whether the parabola -D^(s) curves up or down. The sign of 
{wq — Wp) {w — Wp) determines if the two minima lie at the same side or different sides of the 



^The case Wq — Wp is only possible for A^ < and smoothly connects cases (a) and (b). This case will be 
very relevant in the particular scenario of Section HI The case w — Wp, on the other hand, is only possible 
for Am > and smoothly connects cases (c) and (d). 
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Figure 3: Curves D\{s) (solid blue) and several D'^{s) (red solid, dashed and dash-dotted) 
with different values of ml, intersecting to give two potential minima (indicated by black 
dots) in the 4 different cases listed in eqs. (I2.37p -( l2^40|) . 

pole at s = Wp. Examples for the curves Dl{s) and -D^(s) in these four cases are depicted 
in Fig. [3l Is it obvious that, once we have two degenerate minima in our potential, there is 
a barrier separating them. The two minima and the saddle point in between appear at the 
intersections between the -D?(s) and D'^{s) curves, which must be of sufficiently high degree 
to allow for such structurejj The different cases listed above will in general lead to different 
shapes of the potential barriers, which have an impact on the profile of the critical bubbles 
for the EWPhT. 

It is convenient to introduce a different parametrization for this degenerate case, that can 
be easily connected to the qualitative features of the curves -D^(s) and D'^{s). Remember 
that the solutions of dV/ds = at the axis h = lead to a cubic equation and up to two 
local minima that are cumbersome to determine analytically. To avoid this problem, it is 



^Some analysis in the literature solve the equation dV/ds = for s{h) and then look for a barrier in 
the one-dimensional potential V[h{s), s]. While this is justified in some cases, Fig. [3] illustrates some of the 
possible dangers of this procedure: h{s) might not be single- valued; V[h{s), s] can venture into the unphysical 
region; and sometimes there is no continuous path connecting both minima and having dV/ds = 0. 
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helpful to treat the position of the miniinuni wq as a free parameter and trade it for the 
parameter m*. Imposing the condition that both minima are degenerate, the parameter A^ 
can be fixed. Finally, one can also trade the parameter m^^ for Wp, the point that marks the 
axis of symmetry of the curve D^^ls). We end up with the following parameters 



{w,Wp,Wo,v,ml,ml,Xm} 



related to the ones in eq. (12. 6 p by 



Aw 



m* 



A^ 



mlh 



X2 - 



-X. 

1 



■mt 



+ 



m 



sh 



Aw"^ vAw 



^h 

Aw2 



Aw'' 



mz. 



+ 2^ + 3 



■m 



sh 



vAw 



4 



Am + 



mi 



Aw^ 



H 2 



+ 



ml 



Aw"' 



mi 



Aw'' 



+ 



mt 



+ 2 



m, 



sh 



vAw 



XmV{w - Wp) 



(2.41) 

(2.42) 

(2.43) 
(2.44) 



wherq^ Aw = w—wq and, as promised, we give the expression for A^, the value of A^ required 
for degeneracy of the minima (and already discussed at the end of SubSection l2.2.3p . Some of 
the particular combinations of masses that appear above have a direct physical interpretation. 



For instance, we have 



m%_ 

Aw' 



+ 



y2 



vAw 



Aw' 



m„ 



(2.45) 



where m"^ is the squared mass at the broken minimum along the direction ip = s cos 6 + h sin 6 
(that joins both minima), where the angle 6 satisfies tan^ = v/ Aw. 

Finally, we can now obtain the condition necessary to ensure that wq is the deepest 
minimum along /i = 0, which simply reads: 

,2 



DetA^^ > -% 

Aw^ 



(m^lo)m^ , 



(2.46) 



where rniW^ is the squared-mass d'^V/{dhy at the symmetric minimum (0, Wq), which is given 

by 



^hlo 



Aw^ 



ml 



Am--Pj-2^ >0. (2.47) 

Aw'' vAwj 

It is interesting that these conditions also ensure that the potential is bounded from below: 

that is, if (I2.46P and (I2.47P are satisfied then A^ fulfills the stability constraints discussed in 

Subsection 12X11 

The fact that the quantities in the discussion above depend always on the ratios Rh = 
m,1/Aw'^, Rs = m^/f ^ and Rgh = ^h/i'^^'^) (plus A.^) is a consequence of the fact that, as 



^'^Aw cannot be zero in the presence of a barrier, and therefore, at the EWPhT there is a jump both in 
(h) and (s). This can be important for some EW baryogenesis mechanisms. 
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we mentioned earlier, the conditions for degeneracy are independent of possible rescalings of 
the fields h and s. If such field rescaling is followed by a global rescaling of the potential so as 
to leave A^ unchanged, the above ratios are also invariant under such combination of rescal- 
ings. Multiplying these 4 parameters by a common factor changes the potential by the same 
overall factor and so, the shape of the potential with two degenerate minima is determined 
by just 3 shape parameters: the three independent ratios {\m/ Rh, Rsh/ Rh, Rs/ Rh} ■ 

The potential with such degenerate vacua takes the form 

+ A^slisr - If + A%^isr - IM + sl{2sr - 3)]1 , (2.48) 

where we have used 

, _h _ s — wq 

V Aw 

with the EW minimum at hr = l,Sr = 1 and the symmetric one at hr = Sr = 0. 

The previous discussion assumes Am 7^ 0, but one can also get degenerate minima for 
\m = 0. In that case, the curve -D^(s) becomes a straight line and reads 

2 
h^ = Dl{s) = v^-2v{s-w)-f , (2.50) 



intersecting the axis /i = at 



^h 



.. = »+^. (2.51) 



We can again impose degeneracy of the two vacua, and write m^, in terms of wq, as we did 
for Xm 7^ 0, arriving at the same eqs. (I2.42p and (I2.43p : on the other hand, m^^^ is now fixed 

2{w^ - w) 
The parameters to describe the potential with degenerate minima in the case A^ = are 
then 

{w,W:,„wo,v,ml,ml} . (2.53) 

2.3.1 Flat Directions 

An intriguing situation occurs in the particular limiting case in which D^ ^{s) fall on top of 
each other; then these curves correspond to a flat direction of the potential. In terms of the 
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original potential parameters, -D^(s) = -D^(s) requires 

filfim + ^>^hf^l = , (2.54) 

AXmfil -l^lr + SXhjjl = , (2.55) 

8\hfi3 - ^^mfJ-m = , (2.56) 

A\hXs-\l = 0, (2.57) 

which, in terms of our parameters, simply reaq^ 

A2 = , m, = , DetTW^ = . (2.58) 

While these conditions would be unacceptable at T = 0, we will see in later Sections that, if 
such fiat direction develops at the critical temperature for the electroweak phase transition, 
the strength of this transition can be significantly larger: the small effect from thermally 
induced cubic terms in the finite-temperature potential is enhanced along such fiat directions. 

Imposing the conditions derived above, the potential takes the simple form 

V = -^ [ml{h^ - v^) + X^v\s - wf ± 2mhmsv{s - w)]^ , (2.59) 

where the ± sign corresponds to the possible sign of m1^. The fiat direction will in general 
be parabolic in the [h?, s)-plane, becoming a straight line for Am = 0. We will examine this 
scenario in more detail in the following Sections. 

3 Strong Electroweak Phase Transitions 

The present model can develop very strong phase transitions if the barrier separating the 
broken and symmetric vacua is produced by tree-level effects (as discussed in the previous 
Section) and not by the cubic term resulting from one-loop thermal contributions of bosons 
(which is the most studied mechanism to achieve a first-order phase transition). Indeed, 
the latter are always proportional to the temperature and lead to a critical order parameter 
Vc oc Tf.. In this case the strength of the EWPhT, characterized by fc/Tc, is independent 
of the temperature and typically proportional to small couplings. For a tree-level barrier, 
on the other hand, v^ is proportional to other dimensionful parameters of the potential and 
almost T-independen1o. In this case Vc/T^ can be large for small critical temperatures. 



^^Here there are only three conditions, since the shift symmetry has not been used to fix one of the 
parameters in the original parametrization. In (|2.58p . on the other hand, this degree of freedom disappears 
since the parameters are shift-independent. 

^^Incidentally, due to this property the strength of such transitions will be insensitive to the gauge-fixing 
subtleties discussed in 45_, . 
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3.1 Evolution of Parameters at Finite T 



When the barrier is produced at tree-level, it is enough to include in the one-loop thermal 
potential the leading terms in the high-temperature expansion that lift the minimum in the 
broken phase. In our model, these terms are 






loop 



-Chh^ + -CsS"^ + mgS 



rp2 



(3.1) 



where 



C/i = 77T 



1 

48 
1 



ma 



12 
1 

12 



9g^ + 3/ + 2(6/1? + 12A;, + A^) 
(2Am + 3As) , 

(/is + Aim) • 



(3.2) 



Here g and g' are the SU{2)l and f/(l)y gauge couplings and ht is the top Yukawa coupling. 
Additional particles coupled to the Higgs or the singlet will in general contribute to these 
quantities. At very high temperature the potential is dominated by this contribution, which 
drives (h) — )■ 0, restoring the EW symmetry [46j. On the other hand, the singlet develops a 
thermal tadpole so that (s) — > Sqo = —m^/cs at high TO In the general case, without an 
s — )• — s symmetry, there is no reason to expect s — )■ 0. This could be arranged by using the 
coordinate frame /is = —fim (provided Cg > 0) but there is no sense in which a symmetry 
associated with the singlet is being restored, simply because there is no symmetryl^ 



The key point in our approach is that the terms in (13.11) can be absorbed in the definition 
of T-dependent parameters 



filiT) ^ ^l + c,{T'-T^ 



(3.3) 



rjiz 



where we use a notation in which, when no temperature is indicated for some T-dependent 
quantity, it is implicitly assumed that its value at T^ is meant, e.g. ml = m1{Tc). We can 
then apply the general results on the structure of the potential derived in Section [21 The 
minima are still determined by the curves -D^(s) and -D^(s), which are now T-dependent and 
induce a T-dependence in the location of the minima. We have 



dDl 



dT^ 



Ch_ 



dDl 



dT^ 



4(m3 + Css) 



^■^Under a singlet shift, 7713 — >■ 7713 + c^ct, so that Sqo transforms as it should. 

^*In the case with a Z2-symmetric potential thermal corrections do not break the symmetry and s 
high T, restoring the symmetry in the vacuum. 



(3.4) 



at 
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The curve -D^(s) approaches the axis h = a.s T increases (keeping fixed its symmetry axis 
and without changing its shape); this guarantees that f -^^ at high T. The evolution of 
D'^{s) is more comphcated in generaL We can simphfy somewhat the analysis by choosing 
Urn = through the shift-symmetry, and then we have 



dD^A 



dT^ 



A^ 



_ "^3 \ 
s J 



{^■m = 0) 



(3.5) 



Now, for the cosmological history from Tc to T = to be acceptable, the EW minimum 
must be the global one at T = 0: H(0) < K(0). This requirement will put a constraint on 
the parameters of the potential. The evolution of the difference AVbs(T) = Vb(T) — VgiT) 
with T can be determined as follows: 



dAVbsjT) 
dT^ 



E 



'dV 


dV 


s 



df^ 
dT^ 



+5: 



dV 



d(t)a 



d{(t)a)l 

dT^ 



dV 



dT^ 



(3.6) 



where we symbolically write fi"^ = {/i|,/i^,yuf}, (pa = {h,s}. Noting that dV/dcpa = at 
both minima and using the T dependence of the ^f parameters from eq. (13. 3p . we obtain 



^^^|#^ = \ {CHV\T) + Aw{T) {cAw{T) + w,{T)] + 2m,)] 



(3.7) 



(which is a shift-invariant expression, see footnote [T3|) . A necessary condition for the EWPhT 
to take place is that this derivative is positive at T^ so that the broken minimum is the deepest 
one at least for T ^ Tc, 



dAVUT) 



dT^ 



- {chv"^ + Aw {cs[w + Wo] + 2m3) } > . 



(3.8) 



Note, however, that this is a necessary but not sufficient condition to guarantee that {v,w) 
is the global minimum at T = 0: this must be checked separately, as summarized in Tabled! 
and as will be shown in particular examples in later Sections. 

3.2 Strategy to Identify Strong EWPhTs 



What are the regions of parameter-space that lead to a tree-level barrier? In terms of the 
original parameters of the potential in eq. (12. ip . the answer to this question generally in- 
volves a complicated superposition of non-linear conditions, with hidden physical meaning 
and hard to use for phenomenological applications. This task is greatly simplified by the 
parametrization introduced in Section [21 which allows an easy identification of a potential 
with stable minima. Moreover, when minima exists they must necessarily be separated by 





Parameters 


Conditions 


Degenerate V, (2.48]) 


{w,Wp,WQ, 

'v,ml,m'i,Xm} 


Stability in wq and w, (I2.46p-(l2.47p : 
ml\o,ml,m^^ >0 


Matching 




Broken min. deepest, (13.8]) : 


T<Tc 
General V, (2.14]) 


{v,w,ml,ml, 
"^L)Am,A2,m,} 


V bounded below, (|2.1fip : 

A^ > (A^ < 0) 
A2>-Ay4 (A™>0) 

Vacuum stability. (12.15) : 
DetA^2 > 

Global min., (2.2V]) : 
A2 > 8AV9 



Table 1: Summary on the strategy and parameter conditions to identify potentials with large 
tree-level barriers. 



a barrier. Indeed, any potential of the general form (12.140 . has a stable global minimum at 
{v,w) for any values of the parameters {f , w,m^,m^, m^^, A^, A^, m,,,} if the simple condi- 
tions discussed below eqs. ( 12. ISp . ( 12716]) and ( I2.27P are fulfilled. Similarly, a potential of the 
degenerate form (12.480 . with parameters {w,Wp,wo,v,m'l,m1, Xm} satisfying the conditions 
of eqs. (I2.46p - (l2.47p . has a barrier between two degenerate minima at {v,w) and {0,wo). 
Both cases are summarized in Table [H 

Hence, the strategy to find a model with a strong phase transition is the following, as 
illustrated in Tabled! First, choose a value for the parameters {w,Wp,Wo,v,ml,m'j., Xm}, 
subject to the simple conditions summarized in the upper part of Table [H Any such choice 
determines a potential of the form (I2.48P with two degenerate minima with broken and 
unbroken EW symmetry and a barrier separating them: this will be the thermal potential 
at some critical temperature Tc, which at this point we are free to choose. Once we select 
Tc, we can match this potential to a general potential (12.140 . making sure to satisfy the 
conditions in the middle part of Table [T], which ensure that the broken minimum gets deeper 
than the symmetric minimum for decreasing T < T^. Finally, using the formulae outlined in 
the previous Subsection, we can evolve all the parameters with T to obtain their values at 
T = 0: these are the relevant parameters that enter physically meaningful quantities like the 
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scalar masses, mixings, etc. As T is lowered, it is crucial that our (broken) vacuum remains 
the global stable minimum of the potential: this is guaranteed by the conditions in the lower 
part of Table [H In fact, the stability conditions on A^ are guaranteed to be satisfied once 
the Tc parameters satisfy the conditions in the upper part of Table [H Eventually, h and s 
should be suitably rescaled to ensure that the zero temperature vev is f = vew- In this way, 
different values of Tc will generate a family of models with different values of the potential 
parameters, but all with the same potential shape at T^. 

With the zero temperature potential at hand, a full-fledged one-loop analysis can be 
performed to confirm the first-order nature of the phase transition and to calculate the 
real critical temperature at which the broken and symmetric vacua are degenerate. This 
temperature will in general differ from the T^ parameter we have used, which corresponds 
in the mean-field approximation to the real critical one. We call Tc^mf the latter and Tc^il 
the former. We will show this strategy at work in some examples in later Sections. The 
details of the calculation of the one-loop scalar potential at finite T, which are standard, are 
relegated to Appendix |Al 

Notice that our estimate of Vc/Tc is conservative since the true critical temperature at 
which the transition starts (the nucleation temperature) is smaller than our Tc. For a com- 
plete analysis, one should also recalculate how the sphaleron energy is affected by the Higgs 
barrier and how this impacts the critical ratio Vc/Tc required for a successful preservation 
of the baryon asymmetry. However, it is generically the case that the sphaleron energy is 
dominated by gauge degrees of freedom with Higgs effects amounting to a few percent change 
(see e.g. [S]). 

Before moving to the examples, let us finally mention the case in which the potential 
has two degenerate minima but both with h > 0. Does this correspond to some situation 
of physical interest? We know that, also in such cases, the EW symmetry will be restored 
at some higher temperature and a local minimum at h = will arise, at which point there 
can be only one broken minimum (as we showed in the previous Section). Hence, in this 
case either the two broken minima merge together or one of them moves to h = before 
the critical temperature is reached. Both options correspond to peculiar phase transition 
histories and merit study, which we leave for a future analysis. 
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4 Special Cases: Z2-symmetric Potential 

In a general study like this one, concentrating on the potential (12. ip which involves many 
parameters, it is crucial to identify whether or not some regions of parameter space are more 
natural than others. This point is especially relevant in the presence of symmetries, which 
select a region of parameter space with vanishing volume (and hence unlikely from the point 
of view of a general analysis) and preserve it under RG-evolution and at finite temperature. 
For the SM plus a singlet, the only symmetry (both of the kinetic terms and of the potential 
at renormalizable level) that is interesting from the point of view of the EWPhT is the Z2 
symmetry s — )■ —s. One particular case of interest that falls in this category is the so-called 
Singlet Majoron Model [9]. The EWPhT in this model has been studied in [Hjl^ 

Making the Z2 symmetry manifest (although the general analysis of Section |2] can be 
carried out without problem), the potential is of the form given in eq. (12. ip with 

/il = , flm = , /X3 = . (4.1) 

In terms of our parameters (12. 6p . these constraints translate into two separate branches, 
depending on whether the Z2-symmetry is broken spontaneously or not. The Z2-symmetric 
case has 

w = , m, = , mlf^ = , (4.2) 

and the Z2-broken case has w 7^ and 

V w 
ml = 2\sw'^ = (4A^ + A^) — 2~ ' "^L = ^mVW , m* = X'^w , (4.3) 



ml 



which allows to extract the usual parameters w, m,^, and m^^j in terms of the others: 

w = = , m^ = X = , m,!, = Am — - . 4.4 

V4A2 + XI v^AX' + XI V4A2 + A^ 

This model can then be described by the 5 parameters {f , m\, m^. Am, A^} in both branches 
but, to avoid confusion, the first part of this Section will be clearer in the standard notation 
ofeq. (EH). 

At high temperature we expect the minimum to lie at the symmetric point /i = 0, s = 
but it might happen that s 7^ prior to the EWPhT. The stationary points of the potential 
will be determined by the intersections of the curves dV/dh = and dV/ds = 0, which now 



^^ Although in the model of [Tl] the scalar is complex (it carries lepton number), from the point of view of 
the potential for the real part of s, it reduces to our case. 
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have very simple expressions: 

^ = ^ |/. = 0, and h' = Dlis) = ^^{2^^l - X^s'] 

-^ = ^ js = 0, and h^ = D''^{s) = -^{i4 + XsS^)\ . (4.5) 

OS [ Xm } 

That is, now -D^(s) and -D^(s) are just parabolas with the same axis of symmetry, at s = 0, 
and different width; 



fl 



Furthermore, to understand the nature of stationary points along 



the branch s = 0, it is useful to write 



= ]-X^[h'-hl]. (4.6) 

s=0 ^ 



{ds 

meaning that, for A^ > (A^ < 0) minima along s = can only appear for h"^ above (below) 
the vertex h^ = — 2/i^/Am of the D'^{s) parabola. 

Can this constrained setting give rise to a tree-level barrier? As we saw in the previous 
Section, in order to arrange for two degenerate minima, one at h = and the other at /i 7^ 0, 
no tree- level barrier can appear if the broken minimum has w y^ 0. This is because the Z2- 
symmetry enforces the existence of two minima with s = ±w but this prohibits a minimum 
at h = 0, as discussed below eq. (I2.29p j^1 Hence, the Z2-symmetric case can only have a tree 
level barrier at the critical temperature if the minimum with /i^ > lies at the symmetry 
axis w = Wp = 0. This situation is illustrated by Fig. H] which shows the intersecting curves 
D^ j,(s) in the {h'^/v^, s/wo)-p\aiae (left plot) and the corresponding potential with its barrier 
(right plot). We focus on this particular case in the rest of this Section. 

Following the approach of Section [2l we identify the conditions for a barrier separating 
the broken and unbroken minima. From eq. (I2.29p . we know that along h = 0, stable minima 
(0, ±Wo) require D^^wq) < which leads to 

A^<-%^<0, A^>0, (4.7) 

where Wq, solution of D1{wq) = 0, is given by 

,,,2 /^s rnlJK^v^ - 2mg) . , .^ 

"° = " a: = vHAX^ + xi) • ^'-'^ 

Furthermore, the condition of degeneracy, eq. (I2.43p . imposes 



^^Which parabola is widest depends on the relative size of the two widths Am/(2A/i) vs. 2Xs/Xm and is 
therefore controlled by the sign of A^ . 

""^^In terms of the two parabolas Df^,,{s), it is difficult to arrange that they cut twice, at the broken 
minimum and at the saddle point in between, because they have the same axis of symmetry, at s = 0. 
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Figure 4: Special scenario with Z2 sjTumetry, A^ > and A^ < 0. Left, curves with 
dV/dh = [Dl{s) and h = 0, blue lines] and dV/ds = [D'^{s) and s = 0, red lines] 
intersecting in the minima at {0,±wq) and (^,0), as indicated by the black dots. Right, 
corresponding potential showing the barrier between minima. 



This condition can be rewritten as 



m. 



^u. 



2v A/iA< 



which will be useful later on. 



(4.10) 



In this degenerate case, the potential then takes the form 



1 



2„,2 



V = -rrif^v 



Rh v^w?, 



— + — -1 

V^ Wn 



(4.11) 



t/i '^ <^o \^ <^o 

showing a concrete example in which the overall shape of the potential is controlled by the 
ratio Rg/ Rh- This is now the only relevant shape parameter (as was to be expected starting 
with only 5 d.o.f.s and removing 3 for rescalings, 1 for degeneracy and no shift freedom) and 
it controls the height of the barrier that separates the symmetric and broken minima. For 
comparison with other cases, notice that this degenerate Z2 scenario corresponds in fact to 
the shape parameters Xm/Rh = 1 + IRg/ Rh and Rgh/Rh = 0. 



4.0.1 Case with Flat Directions 



Applying the general discussion of flat directions in Section |2] to the particular case of the 
Z2-symmetric scenario, we see that a fiat direction arises for 

Xh 



Xr 



X — XhXs — T-^m — • 



(4.12) 



If this happens, then the tree-level potential takes the simple form 



V 



\a ik^ + \^/ 



\ ( i,2 I -'- ^rn 2 

:Xh\h + --— s 

^ Ah 



(4.13) 
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Figure 5: Special scenario with Z2 symmetry, A^ > and A^ = showing a flat direction. 
Left, degenerate parabolas Dl^{s). Right, corresponding potential. 



When thinking about further possible symmetries that could enforce such form of the poten- 
tial one should keep in mind that this form is supposed to hold at some critical temperature, 
not at T = 0. While the symmetric form of the quartic couplings would be approximately 
respected by thermal corrections (as quartics have only a logarithmic dependence on T), 
quadratic terms for h and s do evolve differently with temperature and would break that 
symmetry. We will examine this in more detail in the next Subsection. 

There are two qualitatively-different cases depending on the sign of Am- (The limiting 
case Xm = has little interest, as then s and h are completely decoupled from each other.) 
For Am > 0, the flat direction is the parabola 

^ Xm 2 _ ^H^ 
2 Xh Xh 

closed around the origin. Then the potential looks like a Mexican-hat potential, see Fig. |5l 

Of course such potential would not be acceptable at T = (implying in particular a massless 

scalar) but could be of interest at T = T^: the effect of the thermal cubic from gauge bosons 

can be enhanced by the flatness of the potential, leading to a large v{Tc)/Tc. This is confirmed 

by our numerical analysis. As we will see below, when T falls below T^ the minimum (which 

at T = Tc is not located at any precise point along the flat direction) will be driven either 

to /i = 0, s 7^ (a case which does not interest us) or to /i 7^ 0, s = 0, with a big jump in 

v{T,)/T,. 

For Xm < 0, the degenerate parabolas correspond to two flat directions 



h' + 



(4.14) 



1|A, 



f^l 



h^ - -^s^ = ^ , (4.15) 

^ A/i Ah 

running away to infinity, see Fig. [61 The stability of the potential along such directions 
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Figure 6: Special scenario with Z2 symmetry, Am < and A^ = showing two flat directions. 
Left, degenerate parabolas Dlg{s). Right, corresponding potential. 

should be ensured by mass terms or one-loop quartics. This case could be of interest for the 
transition if the fiat directions intersect h = 0, which requires //^ < 0. The thermal evolution 
of the potential in this case will depend crucially on the thermal cubic. 

4.1 Thermal Evolution and EW Phase Transition 



As we saw in Section [31 at leading order in the high-T expansion, only /i^ and /i| are affected, 
according to eqs. fl3.3pJ^i As a result, when T is lowered from Tc down to T = 0, the parabolas 
D^[s) and D^{s) simply drift away from h = keeping their width and symmetry axis fixed. 
As they move at different rates, determined by [see eq. (13. 4p ] 



dDl 



dT^ 






dDl{s) 
dT) 



2c, 



(4.16) 



their relative position can change, together with the location of the minima in the potential. 
This is illustrated by Fig. [TJ which shows snapshots of D^g(s) at different temperatures. 
At very high T the minimum is at (0,0). The Z2 symmetry breaks spontaneously at some 
critical temperature Tz and the two minima at (0, ±wo(r)) move away from the origin as T 
gets lower. Eventually the EW minimum forms and gets degenerate with the Z2-breaking 
minima at T^. For lower temperatures, the EW minimum is the deepest one. 

Let us focus on cases that lead to a tree-level barrier, which, as explained above, require 
w = 0, Am > and A^ < 0. Once we identify the parameters that give such barrier in the 



^® Note that, in the context of the Singlet Majoron Model, right-handed neutrino Yukawas give a significant 
contribution to c^ and the imaginary component of the (complex) singlet contributes to both Cg and Ch, 



increasing the coefficients of As and \m in (|3.2 
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Figure 7: Snapshots of the T-dependent curves -D^(s) (blue lines) and -Df (s) (red lines) 
intersecting at the T-dependent minima (black dots) in a Z2-symmetric scenario with the 
correct cosmological history. The plots are in order of decreasing T, from left to right and 
top to bottom, with T > Tz (upper left); Tz > T > Tc (upper right); T = Tc (lower left); 
and T = Tc/2 (lower right). 



potential, Vc{h, s), of the form (I4.1ip . we still have the freedom to choose Tc and to perform 
the appropriate rescaling to ensure f (0) = vew- To be specific, the potential is 



V = Vc{h, s) - -(Tf - T'){cHh' + Css') 



(4.17) 



As explained in Section 121 we start at T^ with degenerate minima with broken and 
unbroken EW symmetry: Vf,(Tc) = Vs{Tc). As T is lowered we want that the broken minimum 
gets deeper becoming our vacuum, in which case, w{T) will stay at zero for all T < Tc. One 
has 



AVbs{T) = V[viT),0] - V[0,woiT)] 



AX, 



+ 



4A. 



(4.18) 
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so, to end up at the broken minimum at T = we need the conditio: 

Cs V K v^ 

which can be taken as a constraint on the initial parameters v and wq in a specific model 
where Ch and Cg are known constants. Alternatively, we can separate from Cg its As-dependent 
part [see (I3.2p ] as 

Cs = -Xs + 5cs , (4.20) 

where 6cs = A,„/6 here, but in general can include contributions from other particles coupled 
to the singlet. Then, condition (I4.19p translates into a lower limit on A^: 



^s -^ ■^s,min — , 



2cl - XhScs - '2ch\/cl - Xh6cs 
while an upper limit follows from eq. (14. 7p : 



(4.21) 



)2 

Obviously, Xg^min < Xg^max should be satisfied. 

The T-dependence of our parametrization is as follows. In the approximation we work, 
quartic couplings do not depend on the temperature while the rest of parameters do de- 
pend on it, leading to m\{T)^ ^(T) and v{T). This T-dependence can be extracted from 
eqs. (I2.7p -( |2TT3|) . after feeding in them the T-dependent /i^ and /i^. In this way it is straight- 
forward to extract 

v'iT) = v\^ - f T^ , (4.23) 

where we are always implicitly assuming T < T^ (so that the EW minimum is the global 
one) which is the range of interest to run parameters from T^ down to T = 0. From this we 
can already extract the important ratio viT^jTc as 



T, 



c 




Notice, however, that the phase transition cannot be made arbitrarily strong by choosing 
low Tci in which case v{T^ ~ vew-, since then the tunneling probability becomes small 
and the meta-stable symmetric phase will become stable. Moreover, for low T^ the high-T 
approximation breaks down. 



^^In this simple scenario, this condition coincides with the condition derived from d/SVbs{T) / dT'^ = 
\chv'^ — CsWg] /2 > 0, see the discussion around eq. p.7p . 
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Eq. fl4.24p becomes more meaningful when combined with information on the mass spec- 
trum. For the Higgs mass parameter we have 

ml{T) = 2Kv\T) , (4.25) 

where the physical value of the Higgs mass M| = rn\(Q) fixes \h through M| = 2\hv\y^ [here 
we use capital letters for the T = parameters: M^ = m^(0) and M^ = m^(0)]. Similarly, 
the singlet mass is given by 

ml{T) =ml+ Q^^c, - c)j (Tf - T^) . (4.26) 

Notice that the second term is positive for T < Tc due to XhXs < A^/4 and eq. (I4.19p . 
such that the zero temperature mass of the singlet is larger than the one at the critical 
temperature: M^ > ml. From the condition of degenerate minima, eq. (I4.10p . we can obtain 
the additional relation 

m] = ]^v\T,){\^ - 2^hXs) , (4.27) 

arriving at 

So, the singlet mass squared is a simple linear combination of the two mass scales v'^^/ and 
T^ with positive coefficients. For A^ = Xs^min, the coefficient of T^ in eq. 04.28P is zero and Mg 
is independent of T^. In the case A^ = A^/(4A/i) it is the coefficient of v'^^y that cancels, and 
then Ms increases linearly with T^. This is precisely the limiting case with a flat direction 
at Tc (of Am > type) discussed in the previous Subsection. 

To explore in more detail what masses and v(Tc)/Tc are allowed we will proceed as 
follows. Cases with a barrier are simply found by choosing Am > and — A^ in the interval 
(0, A^/4), which in terms of A^ is equivalent to A^ G {0,Xs,max)- The lower part of this 
interval is removed by the condition (I4.2ip . required to guarantee the correct T — > limit, 
leaving only the interval {Xs,min, Xs^max)- Now, for fixed M^ and using Tc as a parameter, 
we can obtain v{T^/Tc and Mg for different choices of Am and A^ in the appropriate ranges 
described above. The results are as shown in Fig. [8], which presents three representative 
cases: a) Mh = 115 GeV with Am = 0.2 and then A, e (0.03,0.09); b) Mh = 115 GeV with 
Am = 1 and A, G (0.09, 2.29); and c) Mh = 200 GeV with Am = 0.5 and A, G (0.09, 0.19). In 
the upper plot, we show the large value of v(Tc)/Tc (which is independent of As) that can be 
obtained as a function of T^. For v{Tc)/Tc > 4 we cannot trust the high-T approximation, 
so we do not explore smaller values of Tc. In the rest of the plots we show the value of 
the singlet scalar mass Mg as a function of the critical temperature Tc for equally spaced 
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Figure 8: Ratio v(Tc)/Tc (upper- left plot) and values of the singlet scalar mass (rest of 
plots) as a function of the critical temperature for the cases a) M^ = 115 GeV, A^ = 0.2; b) 
Mh = 115 GeV, Am = 1; c) M^ = 200 GeV and A^ = 0.5. Different masses correspond to 
different values of A^ G {Xs,min, K,max), with Mg increasing for lower A^. 



values of A^ G {Xs,min, As,mai)- Higher singlet masses correspond to lower values of A^. One 
concludes that very strong first-order EW transitions can be obtained for a wide range of 
scalar masses. 

4.2 One- loop Numerical Analysis 

So far, we have identified choices for the T = parameters that lead to strong electroweak 
phase transitions in the mean-field approximation. It is straightforward to refine these results 
starting from the same tree-level parameters but including in the scalar potential one-loop 
T = corrections and the full one-loop thermal integrals (which correctly take into account 
Boltzmann decoupling effects) further improved by daisy resummation. Details of this stan- 
dard procedure are given in Appendix A. To illustrate the impact of this refinement, we show 
in Fig. [9] the ratio Vc/Tc in the mean-field approximation (blue dashed line) compared with 
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the same quantity calculated with the one- loop thermal potential just described (black solid 
line). The example shown corresponds to Mh = 115 GeV and Ms = 145 GeV. As expected 
from the correct inclusion of Boltzmann decoupling effects, which tend to increase Tc (see 
inset, where Tc^mf corresponds to the critical temperature in the mean- field approximation. 



while Tc^iL ta 
tion is weaker 
baryogenesis. 



ies into account the full one- loop resummed potential), the one- loop transi- 
than in the mean-field approximation but still strong enough to allow for 



We have also examined at one-loop cases with a fiat direction with A^ > 0, case in which 
thermal cubics from bosons play an important role. Although such thermal effects would 
not give rise to strong transitions by themselves (for weak couplings to the Higgs field), in 
the presence of a flat direction the thermal cubic lifts the fiat direction creating a barrier 
between the broken and symmetric minima, and ensures large jump in v{T)/T, leading to 
strong phase transitions at one-loop. Concerning the naturalness of such scenario, notice 
that all that is required to realize it is that quartic couplings satisfy A^ = (which might 
be the result of some symmetry). This makes the parabolas Dl{s) and D'^{s) equally wide 
and, as they are both centered at s = 0, when they shift with temperature, it is guaranteed 
that they will overlap and give a fiat direction at some Tc. 

Our results contradict some claims in the literature concerning this scenario [37], which 
were focused on transitions driven by the thermal cubic. On the other hand, the analysis 
of the Singlet Majoron model in [14| did find strong transitions. One important difference 
between that study and ours is that, for phenomenological reasons (in particular to be able 
to generate nonzero neutrino masses), the scenarios considered in [14)^^ are restricted to 
w ^ 0, which works against the possibility of obtaining a really strong phase transition 
through tree-level barriers. However, loop effects related to sizable Yukawa couplings (to the 
right-handed neutrinos), that we are not discussing, can help in getting strong EWPhTs. 
On the other hand, it is clear that the example in Fig. 9 of [H] corresponds to a one-loop 
deformation of a tree-level case with a nearly fiat direction. One-loop effects shift the EW 
minimum away from w = and lead to a very light scalar, see eq. 14.31 



^°When comparing the two curves in Fig. HI keep in mind that, for a given choice of the model parameters, 
the resulting v{Tc)/Tc at one-loop is displaced to higher Tc and lower v{Tc)/Tc with respect to the tree-level 
value (like the curves' end-points demonstrate). 

^^A direct comparison of results is difficult because the analysis in [T3] is based on scans of the parameter 
space (for instance, from distributions of A^, As and Am we cannot obtain the distribution of A^, which is a 
crucial parameter for EWPhTs) and different mechanisms operate in different regions of parameter space. 
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Figure 9: The ratio Vc/T^ in the mean-field approximation (blue dashed line) and the one- 
loop approximation (black solid line) as a function of the corresponding critical Tc. The 
relation between critical temperatures in the two approximations is shown in the inset. 
Here Mh = 115 GeV, Mg = 145 GeV and Wq = 100 GeV. T^^mf corresponds to the critical 
temperature in the mean-field approximation, while T^^il takes into account the full one-loop 
resummed potential. 

5 Special Cases: A Supersymmetric Example 



A different context in which a particular region of parameter space is selected by some mech- 
anism, is that of supersymmetry. A supersymmetric version of the SM plus singlet model, 
the Next to Minimal Supersymmetric Standard Model (NMSSM), has been considered since 
long ago as suitable to obtain a strong electroweak phase transition with applications to 
electroweak baryogenesis. In particular, ref. [IE] was the first to emphasize in this context 
the relevance of a tree-level cubic term in the scalar potential. 

Here we will examine a particular deformation of this model, the near-to-Minimal Su- 
persymmetric Standard Model (nMSSM) [19], which differs from the NMSSM in having a 
singlet superpotential with a loop-suppressed tadpole and no cubic term. Refs. [16] studied 
the electroweak phase transition in this model finding strongly first-order cases. In the region 
of parameter space examined in [16] the scalar potential reduces to that of a SM Higgs plus 
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a real scalar with a potential of our general form (12. ip but with 



As = , /is = , A„ > 



(5.1) 



It is interesting to note that a shift of the singlet respects the conditions /is = and A^ = 
so that there is freedom to set /i^ = without loss of generality. In terms of the parameters 
of eqs. (12. 7p - (12. 131) these conditions read 



4 



4 



sh 



XmVW 



(5.2) 



and the model can be described by just 5 parameters, 



{v,w,ml,ml,Xm} 



(5.3) 



The model has sufficient structure in the D^ ^{s) functions to allow for tree-level barriers, 



which must necessarily be of type-(d) in the classification of eqs. (I2.37p -( l2^40|) . that is, of those 
with Am > (imposed by supersymmetry) and {wq — Wp){w — Wp) > 0. This last condition 
can be explicitly checked after translating the conditions of eq. (15. ip into the parameters for 
a degenerate minimum (see Table [T]) or, equivalently, imposing the parameter constraints 
needed to have a symmetric minimum at (0, ti^o) degenerate with the broken one at (f,w), 
which now read: 



2„,2 



Wo=W + 



miv 



Imlw 







-1 + W1 + 



ArngW 
nihV 



(5.4) 



and Wp = 0, from which one gets {wq — Wp) {w — Wp) = wwq, which is always positive. Fig. [TOl 
shows an example of such barrier. 



The potential with degenerate vacua takes the simple form 



V 



- 2 2 

-V m^ 



h s — w 

— + 

V^ Wq 



+ |4f-l 



h^{s — w) 

2/.i,2 



V^Wt 







(5.5) 



where we have one single shape parameter: Rg/Rh, as in the Z2-symmetric scenario, but 
now with Xm/Rh = '^Rs/Rh — 1/2 and Rsh/Rh = 1/2. In the particular case ARs/Rh = 1 
the potential (15. 5p develops a flat direction which, in the (/i^/t>^, s/wo)-plane, is a straight 
line through both minima. In terms of the original quartic couplings of the potential this 
scenario corresponds to the extreme limit A^ = 0, that we do not consider further. 
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Figure 10: Scenario with a strong transition in the nMSSM model. Left, intersecting curves 
Dlg{s). Right, corresponding potential 

5.1 Thermal Evolution and EW Phase Transition 

As in previous scenarios, we start from the potential with degenerate minima, which is 
assumed to hold at some Tc, and then use the mean-field approximation to derive the cor- 
responding potential parameters at T = 0. This temperature evolution is given in eq. (13. 3p 
where now the constants Ch, Cg and m^ will also receive contributions from supersymmetric 
particles (if their masses are not much higher than T). As the model can support a strong 
electroweak phase transition without the need of supersymmetric particles coupled sizeably 
to the Higgs (as needed in the MSSM case) we simplify the analysis by assuming that ther- 
mal effects from superpartners are Boltzmann suppressed. If supersymmetric particles do 
not decouple from the thermal plasma their effects can be included and they will only modify 
quantitatively our discussion of strong phase transitions based on the tree-level potential. 

Note that the condition fim = 0, that we have previously imposed using the shift symme- 
try, is respected by thermal corrections (in the mean-field approximation) and implies that 
Dl{s) is centered at Wp = and has the same width at all T's while its vertex h will evolve 



with T as in the Z2 case. Another good property of this choice is that m^ = 0, see (13. 2p . so 
that nl is also T- independent. This makes the evolution of D'^{s) also very simple: 



DM 






(5.6) 



exactly as in the Z2 case, so that also D'^{s) keeps its shape and simply drifts. As usual, 
to have a successful cosmological history, D^^(s) should move in such a way that the EW 
vacuum forms and gets deeper at low temperature. 

The thermal evolution of the potential minimum in this model is more complicated than 
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Figure 11: The ratio Vc/T^ in a nMSSM scenario in the one-loop approximation (black solid 
line) as a function of the corresponding critical T^. The relation between critical temperatures 
in the two approximations is shown in the inset. Here M^ = 135 GeV, Mg = 145 GeV and 
Wo = -400 GeV. 



in the Z2 symmetric case: v{T), w{T) and irtsiT) are determined by the simple equations: 






^l 



mliT)wiT) = f,',{T)woiT) , 



(5.7) 

(5.8) 
(5.9) 



while rn\{T) is trivially related to v{T) by m\{T) = 2\hv'^{T), as usual. This system of 
equations can be solved analytically, although the resulting expressions for f^(T), w{T) and 
ml{T), being solutions of cubic equations, are not very illuminating and we refrain from 
writing them down. 

Once the temperature dependence of all parameters is known, we can relate the param- 
eters at the critical temperature 



to those at zero temperature 



{v,w,ml,ml,Xm} 



{vew, wew, Ml, Mf , \m} 



(5.10) 



(5.11) 
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and, in particular, determine the critical temperature in terms of the physical parameters. 
Some simple relations are 

Mf = m^^^, (5.12) 

wew 

Once we have the T = potential we can perform a refined one-loop analysis including 
thermal decoupling effects as already discussed in previous Sections to confirm the existence 
of strong phase transitions in the regions indicated by the tree-level analysis. An example 
of the results we obtain is given in Fig. [TTl 

6 Special Cases: Light Scalar 

As a final example we consider a realization of the SM plus a singlet with a very light scalar, 
put forward in [41j as a possible way of increasing the strength of the EWPhT. The model 
has 

\rn = , /il = , /i3 = , A^ = , (6.1) 

and fim, the only coupling connecting the s and h sectors, is small: /im = ^mVEw *^ vew- 
The condition yUi = is not respected by thermal corrections and therefore we have to keep 
it nonzero in our discussion, but with /ii(0) = 0. In our parametrization this reads 

A^ = , A^ = , m, = 0. (6.2) 

This case has only five free parameters that we can take as {v, w, ml, m^, fim}- So, this is a 
very constrained scenario and one sees (cf. Tabled]) that the conditions (I2.46l) - (l2.47p cannot 
be fulfilled: it is not possible to have degenerate minima with a barrier in between. This can 
be easily understood by studying the functions Dl^{s), which are straight lines now: 

Dlis) = ^{2f,l-f,^s), (6.3) 

D'M = -4/i±i^, (6.4) 

and, with such simple structure, cannot lead to a tree-level barrier unless a scenario with 
a fiat direction (or close to it) is realized. Such fiat direction requires fimfJ'h = — 4A/j/xf and 
/i^ = 8A/i/x^. To decide how natural this is, we need to examine the thermal evolution of 
these quantities. 
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Figure 12: Snapshots of the T-dependent hnes -D^(s) (blue) and -D^(s) (red) intersecting at 
the T-dependent minimum (black dot) in a scenario with a light singlet having the correct 
cosmological history. The plots are in order of decreasing T, from left to right and top to 
bottom, with T > Tc (upper left); T = Tc (upper right); and T < Tc (bottom). 

6.1 Thermal Evolution and EW Phase Transition 



In this model, fi^ is T- independent in the mean-field approximation (c^ = 0), while fi\ and 
/xf depend on T in the usual way, now with rris = fim/^'^- The evolution of the two curves 
Dl^{s) with temperature is then quite simple: 






Ch_ 



dDt 



dT^ 



1 
3 



(6.5) 



As seen in the [h"^, s)-plane, when T increases, Dl^{s) shift towards the axis h = without 
rotating. In order to have a proper cosmological evolution, starting in the symmetric vacuum 
at high T and ending up at T = in the broken EW minimum, the line -D^(s) should have 
slope smaller than that of D'^{s) (this requires SXhfi'i > /x^) ^^^ it should move faster with 
T^ than D'^{s) (this requires Ch > A/j/3, which is automatically satisfied). An example to 
illustrate this scenario is presented in Fig. [121 In the mean-field approximation the EWPhT 
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Figure 13: The ratio v^/Tc in the scenario with a hght scalar of Section |6] in the one-loop 
approximation as a function of the light scalar mass. Here Mh — Ms2 — 127 GeV and 
Hm = -40 GeV. 



would be second-order. However, a large v{Tc)/Tc could be achieved after including in the 
potential the full one-loop corrections (that include a cubic term for h not only through 
transverse gauge bosons but also through the Higgs itself) if the model is tuned to have a 
nearly flat direction at T^. The tuning involved, which we quantify by the small parameter 
< e ^ 1, requires nearly equal slopes for Z)^^(s): 



SXhul = /^m(l + e) = d(l + e)w|iy • 



In terms of more physical parameters this reads 



4M^M^ ~ etv" 



h s 



m'^EW ' 



(6.6) 



(6.7) 



which indeed requires a light scalar singlet, corresponding to the field excitations along the 
flat direction. Indeed, including the effects of s-h mixing, the two scalar mass eigenvalues 
are 

1 . 7,1... 

(6.8) 



M^c^Ml^ Mlc.\ele'^ 



h 



Although it is difficult to imagine a symmetry reason that could lead to the relation (16. 6p . 
once that tuning is arranged, temperature corrections do not spoil it. In the mean- field 
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approximation the critical temperature is easily computed to be 

which shows that Tc can be made much smaller than in the SM, helping to increase v{Tc)/Tc. 
To determine this quantity one needs to carry out the one-loop analysis. Our renormalization 
conditions for this particular scenario are detailed in Appendix \M Fig. [13] shows our results. 
We have varied e between 0.08 and 0.2, which roughly corresponds to light scalars with 
masses M^^ = 10 -^ 25 GeV, and we have fixed fim = ^40 GeV (this parameter has little 
infiuence on vjT^) and Mh = 120 GeV (leading to M^., ~ 127 GeV). We see that quite 
strong EWPhTs are possible, with Vc/Tc increasing with decreasing e (or M^J, which makes 
the flat-direction flatter. 

7 Numerical Examples. General Case 

From the general analysis of the tree-level potential in Section |2] we have learned that the 
SM with a singlet has a very rich structure. In particular we showed in eqs. (I2.37p - (l2.40p 
that there are four distinct types of arranging for two degenerate minima in the potential, 
one of which, at {v,w), breaks the electroweak symmetry while the other, at {0,Wo), does 
not. The potential with these degenerate vacua is given by eq. (I2.48p . 

7.1 Thermal Evolution and EW Phase Transition 

As already explained in the particular cases studied in previous Sections, we take the previous 
potential with degenerate minima to hold at some critical temperature T^. and use the mean- 
fleld approximation to the free-energy to obtain the corresponding tree-level potential at 
T = 0. The temperature dependence of the potential parameters is very simple but flnding 
how this dependence affects the minima, in particular v{T) and w{T), requires solving a cubic 
equation. A smart choice of the singlet shift can simplify this task by leading to simpler 
analytical expressions for these quantities. Here we simply perform this thermal evolution 
of parameters down to T = numerically. 

We start from one particular potential with degenerate minima, like those just discussed, 
expressed in terms of a set of original parameters {fi\, /z^, fim, /^s, /^i, ^h, ^s, ^m}- If we rescale 
all mass parameters by the appropriate power of some factor A{Tc) we still have one potential 
with degenerate minima corresponding to the set 

{A(r,) VL MTcff^l ^(Te)/i„, A(T,)/X3, ^(Te) V?, K, K, Xm} • (7.1) 
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case 


M,i / GeV 


Ms2 1 GeV 


sin(ahs) 


Am 


w /GeV 


wq 1 GeV 


a) 


115.5 


209.3 


0.612 


-0.1 


600 


120 


b) 


203.9 


114.2 


-0.429 


-0.1 


450 


-90 


c) 


159.2 


108.9 


0.398 


0.2 


100 


-520 


d) 


116.9 


215.5 


-0.621 


0.15 


130 


676 



Table 2: The parameters used in the general models of Fig. [TH 



The thermally corrected potential (in mean-field approximation) for that rescaled potential 
reads 



Vt{Ks) 



^^AiT,)'lxlh' + ^\,h' + ^AiT,)'fj.y + 1a,/ + ^A{n)ixmsh' + ^Xms'h' 



+A{T,ffils + -A{T,)fi3S^ + 



TjChh^ + ^CsS^ + A{Tc)m^s 



(T^ - T^) , (7.2) 



with c/i, Cs and m^ as given in eq. (13. 2p . The factor A{T^ affects all dimensionful parameters 
and is used to guarantee f (0) = vew- The function A{Tc) and the singlet vacuum expectation 
value at T = 0, wewi are obtained by solving numerically the minimization equations 
dVr/dh = and OVt/Os = at T = 0, /i = vew^ s = wew- Once A{Tc) is known as a 
function of T^, eq. (17. 2p describes at T = a family of potentials (parametrized by T^) that 
lead to a strong electroweak phase transition. 

Such potentials can be taken as the starting point for a full one-loop analysis of the 
EWPhT, including one- loop corrections at T = and finite temperature, see Appendix |Al 
We have performed such analysis for some representative families of potentials of the four 
different types discussed in Section |2l For our numerical examples we have chosen the 
parameters as given in Table 121 where Ugh is the scalar mixing angle. The results are 
presented in Fig. [TU which shows the ratio v{Tc)/Tc as a function of Tc both in mean- 
field approximation (dashed line) and one- loop (solid line). The critical temperatures in 
both approximations differ and we also plot their ratio Tc^mf/Tc,il in the insets of each 
plot. Although the generic effect of including the full one- loop thermal corrections is always 
to increase the critical temperature and lower v{Tc)/Tc, Fig. [T3] shows that very strong 
transitions can be obtained, with v{Tc)/Tc > 1 as required for successful baryogenesis, in all 
cases. 
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Figure 14: The ratio Vc/Tc in the mean-field approximation (blue dashed line) and the one- 
loop approximation (black solid line) as a function of the corresponding critical Tc for generic 
scenarios of different types, as indicated. The relation between critical temperatures in the 
two approximations is shown in the insets. 

8 Conclusions 



Our analysis of strong EWPhTs in the SM plus singlet highlights the richness of possibilities 
this simple extension of the SM offers. By relying on a simple mean-field approximation 
to the finite temperature scalar potential and a judicious choice of parametrization we have 
been able to perform a thorough analytical study of strong EWPhTs triggered by tree-level 
dynamics. We have given a strategy, summarized in Table [H to identify the regions of 
parameter space that would lead to such strong EWPhTs. At the same time, our analyt- 
ical approach has improved the understanding of the mechanisms behind such transitions 
allowing us to uncover new scenarios not appreciated before. One interesting example are 
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those transitions that rely on the presence of a flat direction in the potential at the critical 
temperature Tc, mechanism that could operate in many other models besides the SM plus 
singlet one. 

We have computed the important ratio v{Tc)/Tc both in particular realizations of the SM 
plus singlet model previously studied in the literature and in some representative examples 
of the most general model. In models with Z2 symmetry we have determined that a truly 
strong EWPhT based on a tree-level barrier must proceed from an EW symmetric vacuum 
that breaks the Z2 symmetry to an EW broken vacuum that is Z2 symmetric. In other 
particular models and in general, we have shown that v{Tc)/Tc can be easily larger than 1. 
This is a necessary requirement for successful electroweak baryogenesis (switching off in the 
broken phase sphaleron processes that would erase the created baryon asymmetry). This 
jump in v{T) is necessarily associated with a jump in the singlet VEV, which can also be 
relevant in some baryogenesis mechanisms. In addition such strong EWPhTs could also lead 
to a relic stochastic background of gravitational waves. 

We have refined our analysis going beyond the mean-field approximation by including 
one-loop effects, with appropriate renormalization conditions at T = and inclusion of 
thermal effects with daisy resummation and no high-temperature expansions. Although 
v{Tc)/Tc is lowered by such refinement one still finds strong EWPhTs. 

As a byproduct, our parametrization of the scalar potential of the SM plus singlet, which 
allows a good control over its vacuum structure, can be useful also at T = 0. In fact, the 
conditions summarized in the lower part of Table [1] can be applied at T = to guarantee the 
stability of the EW vacuum and then our parametrization can be applied to phenomenological 
analyses allowing a direct control over physical quantities. 
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A Scalar Potential at Finite Temperature 



One-Loop T = Corrections 

The analysis of tliermal corrections to the scalar potential requires for consistency the 
inclusion of one-loop corrections to the potential at T = 0. These are given by the usual 
Coleman-Weinberg correction [50] 



Vi{h,s) 



1 

QAtt'- 



J2NaM'aih,s) 



log- 



Q' 



a. 



(A.l) 



where a runs over all degrees of freedom (counted by Na, which includes a minus sign for 
fermions) with squared- masses M^{h,s) (which depend on the scalar background fields h 
and s), Ca is a constant (equal to 3/2 for scalars and fermions and to 5/6 for gauge bosons) 
and Q is the renormalization scale, that can be fixed e.g. to the top mass. To this we add a 
counter-term potential 

111111 1 

^ 'i ^ t: 'i 'rt O 

(A.2) 
specifically chosen so as to maintain the main properties of the tree-level potentials at T = 
derived in the text. Those potentials have two minima, the one that breaks the electroweak 
symmetry at (vEWyUiEw) and a symmetric minimum at {0,Wo^ew)- The renormalization 
conditions we use cannot be used for potentials without this structure. In order to avoid 
problems with infrared divergent Goldstone contributions, we find convenient to remove the 
Goldstone corrections to Vi in the renormalization conditions that follow (and we indicate 
this by writing Vi). This is simply a change of renormalization conditions and the shift it 
causes in the potential shape is negligible. For alternative treatments of this complication, 
see e.g. [5l]. 



Explicitly, the renormalization conditions that we use are: 



d{Vi + 6V) 



dh 



0, 



d{Vi + 6V) 



ds 







d{Vi + 6V) 



ds 



0, 



(A.3) 



(with the subindex b indicating evaluation at the broken minimum and the subindex s 
evaluation at the symmetric minimum) to ensure that the one-loop minima are still located 
at their tree-level positions; 










dhds 














(A.4) 
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so that the tree-level mass matrix in the broken minimum is not affected at one-loop; then 

d^{Vi + 6V) 



ds^ 







(A.5) 



to have the same singlet cubic coupling; and finally 

{V, + 6V)\b = 0, {V, + 6V)\s = 0, 



(A.6) 



so that the one-loop values of the potential at the minima are the same as the tree-level ones. 
These 9 conditions allow us to determine the 9 counterterms in 6V as (dropping everywhere 
the subindex EW for simplicity) 



5/^h 



2t;2 



dfi] 



[(3 - 5xl)dh, - (1 - xl)dl^ + x^{3x^ - 2)(9^^,J Vi\b 

i\s-4:AVibs^ 



2xi [ 3ds^ - di + -di ]V,\b- 2xi ( ds^V, 



2^^2i(l-3^-) {^dn^-dl)V,\b-^AV,bs +2{l-Qxi)dsy,\ 



+ 



3(1 - 2a:^)(29,, - S^.J + -(1 - 6x^ + 6x^)9^ 



Slim 



v^Aw 
+x^ (2d,^Vi 
3 



^il 



x^ [ 59,, - di + 9d,^ - 2di + -di ) + (1 - 3x^)dl^^ 



Vi\ 



]AV, 



lbs 



(A.7) 



(A.8) 



(A.9) 



2Aw3 
2x^(d,^V,l + 4AVus 



xU5dH^ - dl + 49,J + -(1 - 2x^)at 



Vi\ 



2Aw 



Xq. 



{2-x^)dsMs-K^-Xn,)AV, 



lbs 



(A.IO) 

+ x^{l-x^){bdh^-2dl;)Vi\b 



+ 
1 

^3 



{l-x^ + 2x^)91 + (1 - x^){l - 2x^){di 



L 3 



2a 



Vi\ 



(A.ll) 
(A.12) 



5\s 
5Vo 



v'^Aw'^ 
1 



54. - dl - 3dl,^ 



6ds^ - 2dl + -dt ) \/i|, + 29,,\/i|, - %AV, 



\bs 



2Aw^ 

1 

3 



2 



59,, - 9t - 29t,, + 69., - ^9i 1 \/i|, + 49.,\/i|, - 8A\/i 



16s 



+ 1 -a:^ 



2a:^[(l - xl)d,y^l - 2(2 - x^)A\/i5.] - 4Fi6 
1 



1 



2^3 



:i + a;^)'(54, - 9t) + x^{2 + x^)(29,, - 9^,,) - -x^9: 



(A.13) 
(A.14) 

(A.15) 
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where Aw 



w 



Wo, x^ = w/Aw, AVibs = Vib - Vis and dh^ = vd/dh, ds^ = Awd/ds . 



These counterterms have a finite if — )■ hmit, do not suffer from singularities {Aw cannot 
vanish), and do not spoil the Z2 symmetry when that is a symmetry of the Lagrangian. As 
explained in the text, strong transitions in this Z2-symmetric scenario require w = 0. The 
form of the counterterms is much simpler in that case and can be obtained from the general 
formulas above simply setting to zero x^ and all odd s-derivatives of Vi evaluated at the 
broken minimum. In particular, one gets ^/i^ = ^/^s = ^/^i = 0. 

In the model discussed in Section [6] the tree- level potential does not have two degenerate 
minima but rather a nearly-flat direction, so that we do not use the previous prescription. 
In order to keep the nearly-flat structure at one-loop, we use instead the following renormal- 
ization conditions. For the broken minimum we impose that at one-loop it stays at the same 
tree-level location and with the same spectrum: 

Q Q Q2 Q2 Q2 



dh' ds' {dhy dhds' {dsY 



{Vi + 5V) 







(A.16) 



The (tree-level) line -D^(s) cuts /i = at s = s/i and, to maintain this at one- loop we impose 



diVi + 5V) 



d{h^ 



. 



(A.17) 



(0,Sh) 



Finally, the (tree- level) line D'^{s) cuts h = aX s = and, to maintain this and the same 
slope at one-loop we impose 



d{Vi + 6V) 



ds 



0, 4/.; 



^d\Vi + 5V) 



d{h^)ds 



l^r. 



d^{V^ + 5V) 



ids) 



(A.18) 



where the subindex indicates evaluation at h = 0,s = 0. 



Finite Temperature Potential. 

At the large temperatures of the early Universe plasma, high-temperature effects modify 
the Higgs effective potential (or rather, free-energy). The contribution of the different plasma 
species to the potential, in the non-interacting gas approximation, is given by standard one- 
loop (bosonic/fermionic) thermal integrals. Each particle species, labelled by a, contributes 
to the potential 



6aVT{h, S) 



rp4 poo 



27r2 

T 

12^ 



dx x^ log 







lie 



-^x2+M2(h,s)/T2' 



KbN^{Ml(}i,s) - Ml^^{}i,s,T 



(A. 19) 
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where MT,a{h, s, T) is the thermally corrected mass of the corresponding species and the plus 
(minus) sign in the integrand is for fermions (bosons). The second line in (IA.19I) is present 
only for bosons (we represent this symbolically by writing Sab) and takes into account the 
effect of resumming hard-thermal loops for Matsubara zero modes. For our numerical work 
we used a series expansion of these integrals in terms of modified Bessel functions [26], 
avoiding high-T expansions. 

The thermal masses in the SM plus singlet model are as in the Standard Model except 
for the scalar sector. The squared mass matrix for h and s is 



M' 



-fll + SXhh"^ + i(/im + Xms)s 



■oif^m + 2XmS)h 



l{fi^ + 2Xms)h /i^ + SA^s^ + 2/i3S + iA^/i^ 



1 

+ 48 



9^2 + 3^'2 + 2(6/^2 + 12Xh + X, 






4(2A„ + 3A, 



rp2 



The thermally corrected mass for Goldstones is 



rriG = -f^h + ^hh' + -(/i„^ + Xms)s + — 9^" + ?>g" + 2{Qhi + l2Xh + A 

z 4o 



./2 



T' 



(A.20) 



(A.21) 
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